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Abstract 

The paper begins with a brief recollection of interactions of the 
author with Ya B Zeldovich in the context of the study of relic grav- 
itational waves. The principles and early results on the quantum- 
mechanical generation of cosmological perturbations are then sum- 
marized. The expected amplitudes of relic gravitational waves are 
different in different frequency windows, and therefore the techniques 
and prospects of their detection are distinct. One section of the paper 
describes the present state of efforts in direct detection of relic gravita- 
tional waves. Another section is devoted to indirect detection via the 
anisotropy and polarisation measurements of the cosmic microwave 
background radiation (CMB). It is emphasized throughout the paper 
that the conclusions on the existence and expected amount of relic 
gravitational waves are based on a solid theoretical foundation and 
the best available cosmological observations. I also explain in great 
detail what went wrong with the so-called 'inflationary gravitational 
waves', whose amount is predicted by inflationary theorists to be neg- 
ligibly small, thus depriving them of any observational significance. 



1 Introduction 

The story of relic gravitational waves has revealed the character of Ya. B. Zel- 
dovich not only as a great scientist but also as a great personality. One should 
remember that the beginning of the 1970's was dominated by the belief that 
massless particles, such as photons, neutrinos, gravitons, cannot be created 

1 A contribution to the international conference on cosmology and high-energy astro- 
physics "Zeldovich-90" held in Moscow, 20-24 December, 2004; http://hea.iki.rssi.ru/Z-90/ 
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by the gravitational field of a homogeneous isotropic universe. Zeldovich 
shared this view and was publishing papers supporting this picture. He was 
enthusiastic about cosmological particle creation [1] and contributed a lot 
(together with coauthors) to this subject. However, he thought that some- 
thing interesting and important could only happen if the early universe was 
highly anisotropic. 

When I showed [2, 3] that massless gravitons (gravitational waves) could, 
in fact, be created by the gravitational field of a homogeneous isotropic uni- 
verse, a considerable debate arose around this work. I argued that the cou- 
pling of gravitons to the 'external' gravitational field follows unambiguously 
from the equations of general relativity, and it differs from the coupling of 
other known massless particles to gravity. In contrast to other massless fields, 
this specific coupling of gravitational waves allows their superadiabatic (para- 
metric) amplification by the 'pumping' gravitational field of a nonstationary 
universe. (A similar coupling to gravity can be postulated for the still hy- 
pothetical massless scalar field.) If classical gravitational waves were present 
before the era of amplification, they would have been amplified. But their 
presense is not of necessity: even if the waves are initially in their quantum- 
mechanical vacuum (ground) state, the state will inevitably evolve into a 
multi-particle state. In phenomenological language, gravitational waves are 
being generated from their zero-point quantum oscillations. 

The intense debate has finished in a surprising and very flattering way 
for me. It is common knowledge that it was virtually impossible to win a 
scientific bet against Zeldovich - he knew practically everything about physics 
and had tremendous physical intuition. But sometimes he would find a cute 
way of admitting that his previous thinking was not quite right, and that he 
also learned something from a debate. On this occasion it happened in the 
following manner. 

After one of his rare trips to Eastern Europe (as far as I remember, it was 
Poland) Ya.B. gave me a gift. This was a poster showing a sophisticated, 
impressionist-style, lady. The fact that this was a poster with a sophisticated 
lady was not really surprising - you could expect this from Yakov Borisovich. 
What was surprising and flattering for me was his hand-written note at the 
bottom of the poster. In my translation from the Russian, it said "Thank 
you for your goal in my net". Ya.B. was hinting at my passion for football, 
and he knew that this comparison would be appreciated much better than 
any other. So, this is how a great man admits a clarification of an error; he 
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simply says "thank you for your goal in my net" . 

It was clear from the very beginning of the study of relic gravitational 
waves that the result of amplification of a wave-field should depend on the 
strength and time evolution of the gravitational pump field. We know lit- 
tle about the very early universe these days, even less was known at the 
beginning of the 70's. 

The best thing you can do is to conider plausible models. The simplest 
option is to assume [2] that the cosmological scale factor a(rf) in the expres- 
sion 

ds 2 = a 2 (7?)[-d?7 2 + (Sij + h^dx'dx^ (1) 
consists of pieces of power-law evolution: 

a(v) = (2) 

where l Q and (3 are constants. Then, the perturbed Einstein equations for 
hij(r),x) simplify and can be solved in elementary functions. In particular, 
the intervals of power-law evolution (2) make tractable the effective 'potential 
barrier' a /a in the gravitational wave (g.w.) equation [2]: 



At" + // 



n 2 



a" 



0, (3) 



where ' — d/drj — (a/c)d/dt. 

Using Eq.(2) and the unperturbed Einstein equations one can also find 
the effective equation of state for the 'matter', whatever it is, which drives 
the intervals of a(rj): 

The somewhat strange form of the index 1 + (3 in Eq. (2) was motivated 
by a serious concern of that time - it was necessary to prove that even a 
small deviation from the exceptional law of evolution a(i]) oc rj guarantees 
the effect of g.w. amplification. It is only in this exceptional case that the 
effective potential a /a vanishes, and therefore the superadiabatic coupling of 
gravitational waves to the nonstationary pump field airf) also vanishes. (The 
analogous effective potential is absent in equations for photons, massless 
neutrinos, and some massless scalar particles.) 
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The convenience of the notation utilized in Eq.(2) is that it parameterises 
the exceptional case by (3 = and deviations from this case by a small (3. 
Indeed, it was shown [2] that the amplitude of the generated g.w. mode is 
proportional to small (3; but it is not zero if f3 ^ 0. At the same time, if f3 
is not especially small, the amplitude of the gravitational wave h p (n), soon 
after the beginning of the superadibatic regime and while the wave is still in 
this regime, i.e. before any further processing of the amplitude, evaluates to 




Estimate (5) is approximate (we will be discussing more accurate formulas 
below) but it contains all the necessary physics. The underlying concepts of 
generation and detection of primordial gravitational waves have not changed 
since the first calculations [2, 3], and it is important for our further discussion 
to recall them again. 

To begin with, we note that Eq.(5) is formulated for the dimensionless 
amplitude h of a given g.w. mode characterised by a constant dimensionless 
wavenumber n. (The h(rj) and ji{rj) mode-functions are related by h = 
n/a.) The wavelength A, measured in units of laboratory standards (as 
Zeldovich used to say, measured in centimeters), is related to n by X(rj) = 
2na(rj)/n. It is convenient to use (and we will always do this) such an rj- 
parameterisation of a{jj) that the present-day scale factor is a{jj R ) = 2l H , 
where In = c/H{jj R ) is the present-day value of the Hubble radius. Then, 
uh = 47r is the wavenumber of the waves whose wavelength today is equal 
the prsent-day Hubble radius. Longer waves have smaller n's, and shorter 
waves have larger n's. 

Expression (5) is essentially a consequence of the following two assump- 
tions. First, it is assumed that the mode under consideration has entered the 
superadibatic regime in the past, and is still in this regime. This means that 
the mode's frequency, instead of being much larger than the characteristic 
frequency of the pump field, became comparable with it at some time in the 
past. Or, in cosmological context, the wavelength A (77) of the mode n, instead 
of being much shorter than the instanteneous Hubble radius c/H(n) = a 2 /a', 
became equal to it at some moment of time rji, i.e. Aj = c/ Hi. For the scale 
factors of Eq.(2), this condition leads to {n/n H )\rji\ ps 1. 

Second, we assume that by the beginning of the superadiabatic regime at 
77 = rji, the mode has still been in its vacuum state, rather than, say, in a 
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strongly excited (multi-particle) state. That is, in the language of classical 
physics, the mode's amplitude near r]i was not much larger than hi(n) ~ 



lpi/^i, where Ipi is the Planck length, Ipi = ^jGh/c 3 . This condition on 
the amplitude follows from the requirement that initially there were only the 
zero-point quantum oscillations of the g.w. field, and the initial energy of the 
mode was (1/2)^0^. Because of the condition Aj = c/ H h we can also write 
hi(n) as hi{n) ~ Hi lpi/c. 

The amplitude of the mode, after the mode's entrance to the amplify- 
ing superadiabatic regime, and as long as this regime lasts, remains at the 
constant level hi(n), i.e. h p (n) ~ hi(n). This holds true instead of the adia- 
batic decrease of the amplitude oc l/a(rj) that would be true in the adiabatic 
regime. In general, the quantity Hi is different for different n's: 



Therefore, a specific dependence on n arises in the function hi(n), and this 
is how one arrives at Eq.(5) in a simple qualitative manner. 

Formula (5) gives the evaluation of the primordial (before further pro- 
cessing) g.w. spectrum h p {n). Roughly speaking, the initial vacuum spec- 
trum h v (n) oc n has been transformed into the primordial spectrum h p {n) ~ 
h v (n)n 1+l3 \ where $ characterizes the scale factor of the era when the tran- 
sition from the adiabatic to superadiabatic regime has taken place for the 
given interval of wavenumbers n. However, the same mode n can sooner or 
later leave the amplifying regime and start oscillating again. Obviously, this 
reverse transition from superadiabatic to adiabatic regime is being described 
by the same theory. 

The final amplitudes at some fixed moment of time (for example, today's 
amplitudes) hf(n) are related to the /i p (n)-amplitudes by 



where (3f characterizes the era when the opposite transition from the supera- 
diabatic to adiabatic regime has taken place (this is why the minus sign arises 
in front of 1 + j3f in the exponent). 

The discussed amplitudes h{n) are in fact the root-mean-square (rms) 
amplitudes of the multi-mode field, they determine the mean-square value of 





h f (n) ~ h p {n)n~ {1+p f\ 



5 



the wave field h according to the general formula 

J lb 

It is necessary to say that in the beginning of the 80's, the inflationary 
cosmological scenario governed by a scalar field [4] was gaining popularity. 
Its central element is the interval of deSitter expansion, which corresponds 
to f3 — —2 in Eq.(2) (rj grows from — oo, 1 + /3 < 0) and w = — 1 in Eq.(4). 
By the time of publication of the inflationary scenario, unusual equations of 
state for 'matter' driving the very early Universe, including such exotic ones 
as p = — e, w — — 1, had already been the subject of cosmological research, 
most notably in the work of A. D. Sakharov [5]. 

The g.w. calculations for the special case (3 = —2 were performed in a 
number of papers (see for example [6, 7, 8, 9]). If j3 = —2, the dependence on 
n vanishes in the general Eq.(5), and the primodial (unprocessed) spectrum 
hp{n) becomes 'flat' - that is, n-independent. Ironically, the prospects of 
direct detection of the stochastic g.w. background characterised by the cor- 
responding processed (today's) spectrum had already been explored by that 
time [3]; the processed spectral index for this model is a = 1 in notations of 
that paper. Ref.[3] also suggested the use of cross-correlated data from two 
detectors and touched upon the technique of 'drag-free satellites' that was 
later developed in the Laser Interferometer Space Antenna (LISA). 

The generality of inflationary, quasi-deSitter, solutions was a serious con- 
cern for Zeldovich during long time. He kept wondering about the sensitivity 
of inflationary solutions to the choice of initial conditions. Nobody would 
take the inflationary scenario seriously if it were a very contrived or unsta- 
ble solution. However, it was shown [10] that the inflationary-type evolutions 
are, in fact, attractors in the space of all possible solutions of the correspond- 
ing dynamical system. This decisive property made inflationary evolutions 
more plausible and appealing. 

2 Direct detection of relic gravitational waves 

The spectrum of h rms (v) expected today is depicted in Fig.l (for more details, 
see [11, 12]). Almost everything in this graph is the result of the process- 
ing of the primordial spectrum during the matter-dominated and radiation- 
dominated stages. The postulated 'Zeldovich's epoch' governed by a very stiff 
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effective equation of state, is also present in the graph, as shown by some 
relative increase of power at very high frequencies. The primordial part of 
the spectrum survives only at frequencies below the present-day Hubble fre- 
quency vn pa 2 x 10~ 18 Hz. The available CMB observations determine the 
amplitude and spectral slope of the g.w. spectrum at frequencies around ujj, 
and this defines the spectrum at higher frequencies. 

Figure 1: Envelope of the h rms {y) spectrum for the case (3 = — 1.9 (n = 1.2) 



log tl(V) 




log v, Hz 



The numerical value of h rms at frequencies around is H is determined by 
the numerical value of the observed quadrupole anisotropy of CMB. As will 
be shown in great detail in Sec. 4, it follows from the theory of cosmologi- 
cal perturbations that relic gravitational waves should provide a significant 
fraction of the observed CMB signal at very large angular scales (barring the 
logical possibility that the observed anisotropies have nothing to do at all 
with cosmological perturbations of quantum-mechanical origin). 

In other words, the final theoretical results do not contain any dimension- 
less parameter which could be regulated in such a manner as to make the 
contribution of, say, density perturbations to the quadrupole anisotropy sev- 
eral orders of magnitude larger than the contribution of gravitational waves. 
These contributions must be roughly equal, but the theory cannot exclude 
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that one of them will turn out to be a numerical factor 2-3 larger than an- 
other. Assuming that relic gravitational waves provide a half of the signal, 
one can find from the observed ST/T m 1CT 5 that h rms (u H ) xs 1CT 5 and, 
hence, it follows from Eq.(5) that lpi/l ~ 1CT 5 . 

The slope of the primordial g.w. spectrum is also taken from CMB ob- 
servations. The commonly used spectral index n (we denote it by a Roman 
letter n in order to distinguish from the wavenumber n) is related to the 
parameter (5 appearing in Eq.(5) by the relationship n = 2/3 + 5. The same 
relationship is valid for density perturbations, to be discussed later. The 
current observations [13, 14] give evidence for n pa 1, which corresponds to 
(3 ~ —2. The particular graph in Fig. 1 is plotted for (3 = —1.9, n = 1.2, 
which tallies with the COBE data [15, 16]. (This spectral index n > 1 implies 
that w < — 1, according to Eq.(4). It is not difficult to imagine that such 
an effective equation of state could hold in the very early Universe, if the 
recent supernovae observations hint at the validity of w < —1 even in the 
present-day Universe !) In simple words, the position and orientation of the 
entire piece-wise function h{v) is defined by the known value of the function 
at the point v = uh and the known slope of the function in the vicinity of 
that point. 

Incidentally, the initial quantum vacuum conditions for gravitational waves, 
at all frequencies shown in the graph, are formulated at the 'initial' moments 
of time, when each wavelength of interest was appreciably longer than the 
Planck length. Therefore, the shown results are immune to the short scale 
ambiguities of the so-called 'trans-Planckian' physics (see for example [17]). 
It is a different matter that at some frequencies the initial state is allowed to 
be a somewhat excited state, rather than the pure vacuum state, without run- 
ning into a conflict with the adopted approximation of small perturbations. 
This exotic possibility and the corresponding modifications of the spectrum 
were discussed long time ago [18] (see also a related work [19]). 

The graph in Fig.l shows the piece- wise envelope of the today's spectrum. 
The displayed result is quite approximate. In particular, it completely ignores 
the inevitable oscillations in the spectrum, whose origin goes back to the 
gradual diminishing (squeezing) of quantum-mechanical uncertainties in the 
phases of the emerging waves and the macroscopic manifestation of this effect 
in the form of the standing-wave pattern of the generated field. (This is also 
related to the concept of 'particle pair creation'.) We will discuss these 
spectral oscillations below. 
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Figure 2: S3 LIGO noise curves and the expected sensitivity f2 ~ 10 4 to 
stochastic gravitational waves 




Frequency (Hz) 

Nevertheless, the graph in Fig.l is convenient in that it gives simple 
answers to the most general questions on the amplitudes and spectral slopes 
of relic gravitational waves in various frequency intervals. For example, it 
shows the expected amplitude h rms = 1CT 25 at v — 10 2 Hz. This is the level 
of the signal that we shall be dealing with in experimental programs. In 
terms of the parameter Q gw (u), 

IT 2 / V \ 2 

<W") = (— ) , 

it corresponds to VL gw 10~ 10 at frequency v = 10 2 Hz and in its vicinity. 

Where do we stand now in the attempt of direct detection of relic gravi- 
tational waves? The sensitivity of the presently operating ground-based in- 
terferometers is not good enough to reach the predicted level, but the exper- 
imenters are making a lot of progress. The data from the recently completed 
S3 run of LIGO [20] will probably allow one to reach the astrophysically 
interesting level of Q gw ~ 10~ 4 , as shown in Fig. 2 (courtesy of J. Romano 
and the stochastic backgrounds group of LSC). Fortunately, the projected 
sensitivity of the advanced LIGO (~2011) will be sufficient to reach the re- 
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quired level of h rms ~ 10 25 , Q gw ~ 10 10 , when a month-long stretch of 
cross-correlated data from the two detectors is available. 



Figure 3: Various LISA sources including relic gravitational waves 
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The ESA-NASA space-based mission LISA (~2013) will have a better 
chance to discover relic gravitational waves. Since the expected spectrum has 
larger amplitudes at lower frequencies, the detectability condition is poten- 
tially improving at lower frequencies. In Fig. 3 we show the LISA sensitivity 
in frequency bins of Af = 3 x 10~ 8 Hz, which corresponds to an observation 
time of 1 year. This observation time should make it possible to resolve the 
g.w. lines from thousands of white dwarf binaries in our Galaxy, radiating 
at frequencies larger than 2 x 10~ 3 Hz. By removing the contribution of 
the binaries from the observed records, or by using sophisticated data analy- 
sis techniques without actually removing the contaminating signals from the 
data, one can effectively clean up the window of instrumental sensitivity at 
frequencies above 2 x 10~ 3 Hz. This window in the area of maximal sensi- 
tivity of LISA is shown in the graph together with the expected level of relic 
gravitational waves in that window. 
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3 Indirect detection of relic gravitational waves 
via CMB anisotropics and polarisation 

The expected amplitudes of relic gravitational waves reach their highest level 
in the frequency interval of 1CT 18 — 10~ 16 Hz. This is why one has very 
good prospects for indirect detection of relic gravitational waves through the 
measurements of anisotropies in the distribution over the sky of the CMB 
temperature and polarisation. (For an introduction to the theoretical tools 
of CMB physics, see for example [21].) 

The accurately calculated power spectrum h^ ms (n) is shown in Fig. 5 [22]. 
The spectrum is calculated at the moment of decoupling (recombination) of 
the CMB, with the redshift of decoupling at z^, c = 1100. The derivation of 
the spectrum takes into account the quantum-mechanical squeezing of the 
waves' phases, which manifests itself macroscopically in the standing-wave 
character of the generated gravitational waves. From the viewpoint of the 
underlying physics, it is this inevitable quantum-mechanical squeezing that 
is responsible for the oscillations in the power spectrum. 

The displayed spectrum was obtained under the assumption that f3 = 
—2 (n = 1), i.e. for a flat primordial spectrum. The survived part of the 
primordial flat spectrum is seen on the graph as a horizontal part of the 
curve in the region of very small wavenumbers n. The normalisation of the 
spectrum is chosen in such a way that the induced quadrupole anisotropy of 
the CMB today is at the level of the actually observed quadrupole [15, 13]. 
Specifically, the temperature function 1(1 + 1)Q in Fig. 4, calculated from 
the spectrum in Fig. 5, gives the required value of 960 (fJ.K) 2 at I = 2. The 
distribution of other induced multipoles is also shown in Fig. 4. 

Figures 4 and 5 are placed one under another on purpose. This place- 
ment gives a better visual description of the fact noticed and explained previ- 
ously [23]. Namely, the oscillations in the metric (gravitational field) power 
spectrum are entirely responsible for the oscillations in the angular power 
spectrum of the CMB temperature, with almost universal correspondence 
between extrema in the wavenumber space n and extrema in the multipole 
moment space I. If there is much/little power in the gravitational field per- 
turbations of a given interval of wavelengths, one should expect much/little 
power in the temperature fluctuations at the corresponding angular scale. 

It is the oscillations in the metric power spectrum that are responsible 
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for the oscillations in the /-space, and not the mysterious explanations of- 
ten repeated in the literature, which claim that the peaks in the function 
/(/ + l)Ci arise because of some waves being caught (at the moment of de- 
coupling) in their maxima or minima, while others are not. To illustrate 
the role of standing gravitational waves and the associated power spectrum 
oscillations, versus traveling gravitational waves with no power spectrum os- 
cillations, it was explicitely shown [23] that the later hypothesis does not 
produce oscillations in the /-space. 

Incidentally, it was argued [23] that in the case of density perturbations, 
the main contribution to the peaks in the temperature function /(/ + 1)C/ can 
also be provided by oscillations in the metric power spectrum, rather than 
by the temperature variations accompanying sound waves in the photon- 
electron-baryon plasma at the last scattering surface. In the case of density 
perturbations, the metric power spectrum is mostly associated with the grav- 
itational field of the dark matter, which dominates other matter components 
in terms of the gravitational field. 

Oscillations in the metric power spectrum in the early Universe are in- 
evitable, and for the same reason as in the g.w. case, namely, because of the 
standing-wave pattern of the metric perturbations that is related to their 
quantum-mechanical origin. Therefore, the often-discussed "acoustic" peaks 
in the /-space may well turn out to be "gravitational" peaks. It remains 
to be seen how this circumstance can change inferences about cosmological 
parameters. 

We shall now turn to the CMB polarisation. (For some important pa- 
pers on CMB polarisation, see for example [24, 25, 26, 27, 28, 29, 30].) It 
follows from the radiation transfer equations that the polarisation of CMB 
is mainly determined by the first time-derivative of the metric perturbations 
in the interval of time when the polarisation is mainly produced. Therefore, 
it is the power spectrum of the function ^-(77, x) that is of a primary im- 
portance. Since the g.w. field itself, including its normalisation, has been 
fully determined, the quantity of our interest is directly calculable. In Fig. 
7 we show [22] the power spectrum (h' / n)^. ms (n) , calculated at the time of 
decoupling. The induced E and B components of polarisation are shown 
in Fig. 6. This graph was derived under the usual assumptions about the 
recombination history, which means, in particular, that the polarisation was 
primarily accumulated during a relatively short interval of time around z& ec . 

Similarly to the case of temperature anisotropies, the extrema in the 
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Figure 4: CMB temperature angular power spectrum 
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gure 5: Power spectrum of gravitational waves at decoupling 




wavenumber, n 
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Figure 6: Angular power spectrum for CMB polarisation 
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Figure 7: Power spectrum of first derivative of the g.w. metric 
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graphs of Fig. 7 and Fig. 6 are in a good correspondence with each other. If 
there is not much power in the first time-derivative of the metric, you should 
not expect much power in the polarisation at the corresponding angular scale. 
On the other hand, the region of wavenumbers n ~ 90, where there is the 
first pronounced peak in Fig. 7, is fully responsible for the first pronounced 
peak in Fig. 6 at the corresponding angular scales I ~ 90. 

Figure 8: Expected numerical level of anisotropy and polarisation induced 
by relic gravitational waves 



10 4 
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In Fig. 8 we combine together some of the expected signals from relic 
gravitational waves. They are encoded in the CMB anisotropies and po- 
larisation. This figure includes also a possible polarisation bump, discussed 
previously by other authors, at very small Vs. This feature arises because 
of the extended reionisation period in the relatively late universe, around 
z re i ~ 17. In agreement with the explanations given above, the amplitude 
and position of this bump in the /-space are determined by the amplitude and 
position of the first maximum in the power spectrum (h'/n) 2 of the function 
h'i^rj, x), calculated at z rei . 

The resulting graphs in Fig. 6 and Fig. 8 are qualitatively similar to 
the graphs derived by other authors before us. However, we take the re- 
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sponsibility of claiming that the numerical level of, say, the B component of 
polarisation shown in our graphs is what the observers should expect to see 
on the sky. Of course, this statement assumes that the observed large-scale 
anisotropies of CMB are caused by cosmological perturbations of quantum- 
mechanical origin, and not by something else. 

The true level of the B signal can be somewhat higher or somewhat lower 
than the theoretical level shown in our figures. But the signal cannot be, 
say, several orders of magnitude lower than the one shown on our graphs. In 
contrast, the inflationary literature claims that the amount of "inflationary 
gravitational waves" vanishes in the limit of the flat primordial spectrum 
/3 = —2 (n = 1). Therefore, the most likely level of the B mode signal 
produced by "inflationary gravitational waves" is close to zero. This would 
make the detection impossible in any foreseeable future. It is a pity that 
many of our experimentall colleagues, being guided by the wrong theory, are 
accepting their defeat even before having started to build instruments aimed 
at detecting relic gravitational waves via the B component of polarisation. 

Their logic seems to be the following: 'we would love to discover the 
fundamentally important relic gravitational waves, but we were told by in- 
flationists many times that this is very unlikely to happen, so we agreed to 
feel satisfied even if we succeed only in putting some limits on, say, polar- 
isation properties of dust in the surrounding cosmos'. The author of this 
contribution fears that in a complex experiment like the B-mode detection, 
this kind of logic can only lead to overlooking the important signal that the 
experiment originally targeted. 

Concluding this section I would like to say as a witness that Zeldovich 
suggested using the CMB polarisation as a g.w. discriminator, as early as in 
the very beginning of the 80's. This was clearly stated in private conversa- 
tions, but I am not aware of any written records. 

4 The false "standard inflationary result" . How 
to correctly quantise a cosmological har- 
monic oscillator 

Why bother about relic gravitational waves if inflationists claim that the 
amount of relic gravitational waves (inflationists and followers call them "in- 
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flationary gravitational waves") should be zero or almost zero? This claim 
is a direct consequence of the so-called "standard inflationary result" , which 
is the main contribution of inflationary theorists to the subject of practical, 
rather than imaginary, cosmology. 

In the inflationary scenario, the 'initial' era of the universe expansion is 
driven by a scalar field with the scalar field potential V(tp). It is in this era 
that the initial quantum vacuum conditions for cosmological pertubations 
are being formulated. The inflationary solutions for the scale factor a{rj) 
are close to the deSitter evolution characterised by (3 = —2 in Eq.(2). The 
effective equation of state for the scalar field is always e + p > 0, so that for 
the power-law intervals of expansion driven by the scalar field, the parameter 
(3 can only be (3 < —2, see Eq.(4). Therefore, one expects the primordial 
spectrum of the generated metric perturbations to be almost flat, i.e. the 
primordial spectral index n should be close to n = 1, with n < 1. 

The beginning of the amplifying superadiabatic regime for the given mode 
of perturbations is often called the 'first Hubble radius crossing', while the 
end of this regime for the given mode is often called the 'second Hubble 
radius crossing'. The "standard inflationary result" is formulated for cosmo- 
logical perturbations called density perturbations (scalar, S, perturbations) 
as opposed to the gravitational waves (tensor, T, perturbations) considered 
in Sec. 1. 

The "standard inflationary result" states that the final (second crossing, 
/) amplitudes of quantities describing density perturbations are related to 
the initial (first crossing, i) values of if and other quantities, according to the 
evaluation: 

(7)r (W '~ (0 '~M?H£)r^ (6) 

The numerator of the last term on the r.h.s. of Eq.(6) is the value of 
the Hubble parameter taken at the moment of time when the given mode 
enters the superadiabatic regime. This is the same quantity which defines 
the g.w. ('tensor') metric amplitude, as described in Sec. 1. Since we are 
supposed to start with the initial vacuum quantum state for all cosmological 
perturbations, one would expect that the results for density perturbations 
should be similar to the results for gravitational waves. One would expect 
that the amplitude hs of the generated 'scalar' metric perturbations should 
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be finite and small, and of the same order of magnitude as the amplitude hr 
of 'tensor' metric perturbations. 

However, according to the "standard inflationary result" , this is very far 
from being the case. The denominator of the last term of Eq.(6) contains a 
new factor: \f\ — n. This factor goes to zero in the limit of the most interest- 
ing and observationally preferred possibility of the flat (Harrison-Zeldovich- 
Peebles) primordial spectrum n = 1. Correspondingly, the amplitudes of 
the generated density perturbations go to infinity, according to the predic- 
tion of inflationary theorists, in the limit of the flat spectrum. (By now, the 
"standard inflationary result" (6) has been cited, used, praised, reformulated, 
popularised, etc. in hundreds of inflationary publications, so it has become 
'accepted by way of repetition'.) 

As will be demonstrated below, the divergence in Eq. (6) is not a violation, 
suddenly descending upon us from the 'blue sky', of the adopted approxima- 
tion of small linear perturbations. This is a manifestation of the incorrect 
theory. Even if the spectral index n is not very close to 1, and you combine 
n with a reasonable Hi in order to obtain, for example, a small number 10 -5 
for the r.h.s. of Eq.(6), this will not make your theory correct. This will be 
just an acceptable number accidentally following from the wrong formula. 
You will have to pay a heavy price in some other places. 

An attempt to derive physical conclusions from this formula can only lead 
to mistakes. The current literature is full of incorrect far-reaching physical 
conclusions derived from this wrong theory. This is a kind of situation which 
L. D. Landau used to describe sarcastically in the following words: "If you 
assume that the derivative of the function sin a; is In re, rather than cosx, you 
can make many wonderful discoveries....". 

In inflationary literature, the 'zero in the denominator' factor \/l — n 
appears in many different dresses. It is often written in equivalent forms, 
such as (ip/H) v (V ip /V) i , {H y(p /H) v ^1 + u>j, etc. Inflationists are routinely 
hiding their absurd prediction of infinitely large amplitudes of density per- 
turbations that should take place in the limit of the flat spectrum n — > 1. 
They divide the g.w. amplitude hr over the predicted divergent amplitude 
hs- This division produces the so-called 'tensor-to-scalar ratio', or 'consis- 
tency relation': hx/hs ~ y/l — n. The quantity common for the T and 
S perturbations, cancels out in the composed ratio, and the 'zero in the de- 
nominator' factor is transferred to the numerator of the final expression. It 
is then declared that the metric amplitude hs of density perturbations is de- 
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termined by the observed CMB anisotropies, and, therefore, the inflationary 
'consistency relation' demands that the g.w. amplitude hx must vanish in 
the limit of n — > 1. 

In other words, instead of being horrified by the fact that their theory 
predicts arbitrarily large amplitudes of density perurbations (and, hence, the 
theory is in complete disagreement with observations, because the analysis 
of the data shows no catastrophic increase in the amplitude when the tested 
spectral index approaches n = 1), supporters of the inflationary approach to 
science systematically claim that their theory is in 'spectacular agreement' 
with observations, and it is gravitational waves that should vanish. 

If this were true, there would not be much sense in attempting to detect 
primordial gravitational waves, as the observations persistently point toward 
n PS 1, including n = 1. It is quite common to hear these days enthusiastic 
promises of inflationary believers to detect "inflationary gravitational waves" 
in the "not-so-distant future" via the measurement of 5-mode polarisation 
of CMB. But from other papers of the same authors it follows that there is no 
reason even to try. If you trust and cite inflationary formulas, the expected 
amount of "inflationary gravitational waves" should be very small or zero. 
You can only hope to be extremely lucky if you suggest to detect them, even 
in the quite distant future, for example with the proposed mission called Big 
Bang Observer. And nobody should be surprised if you have found nothing, 
because n = 1 is in the heart of all claims, theoretical and observational. 
Moreover, most loyal inflationists would say that this was exactly what they 
had always been predicting. 

To demonstrate the incorrectness of inflationary conclusions, we shall now 
concentrate on the 'zero in the denominator' factor. We will have to recall the 
quantisation procedure for gravitational waves and density perturbations. It 
is necessary to remind the reader that some inflationists and their supporters 
insisted for many years on the claim that the dramatic difference in the final 
numerical values of and h $ arises not because of the initial conditions, 
but because of the subsequent evolution. 

Specifically, they claimed that the classical long-wavelength 'scalar' met- 
ric perturbations are capable of experiencing, in contrast to gravitational 
waves, a "big amplification during reheating" . (For a critical discussion, see 
[31].) But it now looks as if the fallacy of this proposition has become clear 
even to its most ardent proponents. Therefore, we shall now focus on the 
issue of quantum mechanics and initial conditions. 
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The perturbed gravitational field for all three sorts of cosmological per- 
turbations (scalar, vector, tensor) is described by Eq.(l). For simplicity, we 
are considering spatially flat cosmologies, whose spatial curvature radius is 
infinite. However, if the spatial curvature radius is finite but, say, only a 
factor of 10 longer than Ijj, very little will change in our analysis. 

The metric perturbations hij(r),x) can be expanded over spatial Fourier 
harmonics labeled by the wavevector n: 



hij(ri,x)= (7) 
; ^, /n / d 3 n V JL(n)— L= 



The three sorts of cosmological perturbations are different in that they have 

s 

three different sorts of polarisation tensors Pij{n), and each of them has 
two different polarisation states s — 1,2. The 'scalar' and 'vector' metric 
perturbations are always accompanied by perturbations in density and/or 
velocity of matter. The normalisation constant C is determined by quantum 
mechanics, and the derivation of its value is one of the aims of our discussion. 

Let us recall the procedure of quantisation of gravitational waves. Let 
us consider an individual g.w. mode n. The time-dependent mode functions 

s 

hn(v) can be written as 

L(V) = -^tUv)- (8) 
a(ri) 

For each s and n, the g.w. mode functions /j,(rj) satisfy the familiar equa- 
tion (3). 

The action for each mode has the form 

S = j Ldr), (9) 

where the g.w. Lagrangian L is given by the expression [32] 

(10) 



L gw - 7^n-*a 2 



■ — n ' 



and 

C 4 
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The Euler-Lagrange equation 

dL d dL 
dh drj dh' 



= 



for the dimensionless g.w. variable h 
motion 

h" + 2-ti - 

a 



fi/a brings us to the equation of 



n 2 h = 0, 



(11) 



which is equivalent to Eq.(3). 

In order to move from 3-dimensional Fourier components to the usual 
description in terms of an individual oscillator with frequency n, we will be 
working with the quantity h introduced according to the definition 



h = 



a 



h 



\fha§ 



-h 



I h Aq 



(12) 



where a is a constant. This constant a is the value of the scale factor 
a{rj) at some instant of time r] = i] where the initial conditions are being 
formulated, and Ao = lixa^jn. 

In terms of h, the Lagrangian (10) takes the form 

2 



1 

2n \a 



212 

n h 



(13) 



The quantity h = q is the 'position' variable, while the canonically con- 
jugate 'momentum' variable p is 



p = ?£. = l(±) 2 h' 

dh' n \a J 



(14) 



In the distant past, at times near 770, and before rji when a given mode 
entered superadiabatic regime, the g.w. amplitude behaved according to the 
law 

1 



71(77) (X 



a(rj)) 



-in{i)— rjo) 



The time- derivative of a(i]) can be neglected, i.e. a' /a n. Then, we 
promote q and p to the status of quantum-mechanical operators, denote 
them by bold-face letters, and write down their asymptotic expressions: 



% ao 
2~a 



ce -m(??-?7o) _|_ c t e »"(>7-»7o) 



(15) 
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_ ce -in( V -V0) + c t e in( V ~V0)] _ 

The commutation relationships for the q, p operators, and for the anni- 
hilation and creation operators c, c^, are given by 



c, c" 1 " 



= 1. 



[q, p] = ih, 

The initial vacuum state |0) is defined by the condition 

c|0) = 0. 

This is indeed a genuine vacuum state of a simple harmonic oscillator, which 
gives at r\ = r] the following relationships 

(0|q 2 |0) = (0| P 2 |0) = 5 AqAp=| 



The root-mean-square value of q in the vacuum state is q rms = yh/2. 
Combining this number with the definition (12) we derive 

*™ = «o|h'|o» 1/a = ^Mpa. (17) 

Extrapolating the initial time t]q up to the boundary between the adiabatic 
and superadiabatic regimes at r\ = rji, we derive the estimate h rms ~ lpi/\- 
It is this evaluation that was used in [2] and in Sec. 1. More accurate 
calculations along these lines produce C = y/TEnlpi in expression (7) for 
gravitational waves. 

A consistent formal derivation of the total Hamiltonian, including the 
terms describing interaction of the oscillator with external field, is presented 
in Ref.[34] by equations (19)-(24) there. Technically, the derivation is based 
on the canonical pair q = /i, p = dL/d/i'. The Hamiltonian associated with 
the Lagrangian (13) has the form 

Hfa) = nc f c + trc t2 + a*c 2 , (18) 

where the coupling to the external field is given by the function 17(77) = 
(i/2)(a'/a). In the same Ref.[34] one can also find the Heisenberg equations 
of motion for the Heisenberg operators c(r/), c'(?7), and their connection 
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to classical equation (3). The asymptotic expressions for the Heisenberg 
operators, 

c(r]) = ce- m ^°\ c^rj) = c ^e in ^°\ 

enter into formulas (15), (16). Clearly, the vacuum state |0), defined as 
c(?7)|0) = 0, minimises the oscillator's energy (18). 

A rigorous quantum-mechanical Schrodinger evolution of the initial vac- 
uum state of cosmological perturbations transforms this state into a strongly 
squeezed (multi-particle) vacuum state [32], but we focus here only on the 
initial quantum state, which defines the quantum-mechanical normalisation 
of our classical mode functions. 



We shall now switch to density perturbations. 



For each mode n of density perturbations (S-perturbations), the mode's 
metric components hij entering Eq.(l) can be written as 

Kj = h(7])QSij + hi(r])n~ 2 Q :ij , 

where the spatial eigen-functions Q are Q = e ±mx . Therefore, the metric 
components associated with density perturbations are characterised by two 
polarisation amplitudes: h(i]) and hi(i]). If the initial era is driven by an arbi- 
trary scalar field ip, there appears a third unknown function - the amplitude 
(fii(rj) of the scalar field perturbation: 

^ = ip (v) + tpi(v)Q- 

One often considers the so-called minimally coupled to gravity scalar field 
tp, with the energy-momentum tensor 



Tp, = WW - 5V 



^ a/3 ww + V(p) 



The coupling of scalar fields to gravity is still a matter of ambiguity, and the 
very possibility of quantum-mechanical generation of density perturbations 
relies on an extra hypothesis, but we suppose that we were lucky and the 
coupling was such as we need. The three unknown functions h(i]), hi(r]), 
Pi(i]) should be found from the perturbed Einstein equations augmented by 
the appropriate initial conditions dictated by quantum mechanics. 
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It is important to note that inflationary theorists are still struggling with 
the basic equations for density perturbations. In inflationary papers, you 
will often see equations containing complicated combinations of metric per- 
turbations mixed up with the unperturbed and/or perturbed functions of the 
scalar field tp and V((p). 

Inflationists are still engaged in endless discussions on the shape of the 
scalar field potential V((p) and what it could mean for countless inflation- 
ary models. However, this state of affairs is simply a reflection of the fact 
that the equations have not been properly transformed and simplified. Since 
the underlying physics is the interaction of a cosmological harmonic oscilla- 
tor with the gravitational pump field, mathematically the equations should 
reveal this themselves. And indeed they do. 

It was shown in paper [33] that, for any potential V(ip), there exists only 
one second-order differential equation to be solved: 



// + /! 



n 



(«Vt)" 



o, 



(19) 



where the function \x{r\) represents the single dynamical degree of freedom, 
describing S-perturbations. The effective potential barrier (a^/7)" '/(a^/7) de- 
pends only on a{r\) and its derivatives, in full analogy with the g.w. oscillator, 
Eq.(3). The time-dependent function 7 (7(77) or 7(t)) is defined by 



H 2 



As soon as the appropriate solution for ^(77) is found, all three functions 
describing density perturbations are easily calculable: 



h(r)) = -H(rj) 



rv 

/ Hyftdr) + Ci 

J 1)0 



h " - ^-h' + n 2 h 



v/7 



-h 



OV7 



The constant Cj reflects the remaining coordinate freedom within the class 
of synchronous coordinate systems. (Another constant comes out from the 
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integration of the above-given equation for hi'.) The funcion \x does not 
depend on this remaining coordinate freedom, and the constant Cj cancels 
out in the expression defining \i{r\) in terms of h{r\): 

" h-^h'. 



H> 

The function is that part of the scalar metric amplitude h(rj) which 

does not depend on the remaining coordinate freedom ('gauge-invariant' met- 
ric perturbation). 

In the short- wavelength regime, the function \i describing density per- 
turbations behaves as ji oc e~ mv . This is the same behaviour as in the case 
of the function fi describing gravitational waves. This similarity between 
the respective functions /i (/i-p and fis) is valid only in the sense of their 
asymptotic 77-time dependence, but not in the sense of their overall numeri- 
cal normalisation (see below). 

In the long-wavelength regime, the dominant solution to Eq.(19) is /i oc 
a^. The quantity which remains constant in this regime is It is 

this physically relevant variable that takes over from the analogous variable 
h = fi/a in the g.w. problem. We introduce the notation 

" C, (20) 



where fi satisfies Eq.(19). 

To make contact with earlier work, it should be mentioned that the pre- 
viously introduced quantity 

Cbst JW^it + #i 

3 1 + w 

where $ is Bardeen's potential and BST stands for Bardeen, Steinhardt, 
Turner [35], can be reduced to our variable ( (20) up to the numerical coeffi- 
cient — (1/2). Our quantity fi for density perturbations can also be related to 
the variable uclms, where CLMS stands for Chibisov, Lukash, Mukhanov, 
Sasaki [36, 37, 38]. 

In preparation for quantisation, we should first identify the inflationary 
'zero in the denominator' factor. The unperturbed Einstein equations for the 
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coupled system of gravitational and scalar fields require [33] 




Therefore, the 'zero in the denominator' factor 

is expressed in the form of very small values of the dimensionless function 
a/7- 

Within the approximation of power-law scale factors (2), the function 7 
reduces to a set of constants, 

2 + 13 2 

7= tt^ 1+w =r- 

The constant 7 degenerates to zero in the limit of the evolution law with (5 = 
—2; that is, in the limit of the gravitational pump field which is responsible for 
the generation of primordial cosmological perturbations with flat spectrum 
n = 1. So, we are especially interested in the very small values of y/j. 

It was shown [33] that the dynamical problem for the scalar-field-driven 
S-perturbations can be obtained from the dynamical problem for gravita- 
tional waves by simple substitutions: a(i]) — > 0(77)^/7(77), I^t(v) ~^ Hsiji)- 
(This is not a conjecture, but this is a rule whose validity was established 
after a thorough analysis of these two problems separately.) Each of these 
substitutions is valid up to an arbitrary constant factor. Using these substi- 
tutions, one obtains the S-equation (19) from the T-equation (3), and one 
obtains the physically relevant variable ( = /Js/a^/7 for S-perturbations from 
the g.w. variable h — Ht/cl- 

Moving from the 3-dimensional Fourier components of the field ( to an 
individual oscillator with frequency n, we introduce the quantity ( according 
to the same rule (12) that was used when we introduced h. Namely, we 
introduce 

7 _ Q o a _ v^QQ a / jj A o a (21) 
^ V^T y/tol Pl n V 327r 3 /p^' 1 ' 
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The application of the substitutions a — > a = h — > £ to the g.w. 

Lagrangian (13) gives rise to the Lagrangian Ld p for the single dynamical 
degree of freedom describing S-perturbations: 

L* = i(^) 2 [(c') 2 -« 2 c 2 

Obviously, the Euler-Lagrange equation 

C" + 2 { ^H, £ + n 2 ( = (23) 

derivable from the lagrangian (22) in terms of the independent variable (, is 
equivalent to Eq.(19) which is the Euler-Lagrange equation derivable from 
the Lagrangian (22) in terms of the independent variable ps- The Lagrangian 
(22) should be used for quantisation. The Lagrangian itself, as well as the 
action and the Hamiltonian, does not degenerate in the limit 7 — > 0, i.e., in 
the limit of the most interesting background gravitational field in the form 
of the de-Sitter metric, 7 = 0. 

We shall start with the analysis of the paper [39] which, together with 
Ref. [40], is sometimes referred to as the most recent work that contains a 
rigorous mathematical derivation of the "standard inflationary result" . The 
author of these papers uses slightly different notations, such as a 2 = e 2p and 
ip — (f). In his notation, the quantity (fio/H is <p/p, so that the 'zero in the 
denominator' factor appears as 0*/p*, where the asterisk means "the time of 
horizon crossing" . 

As a "useful example to keep in mind" for quantisation of density pertur- 
bations, the author suggests the artificial model of a test massless scalar field 
/ in deSitter space. But the Lagrangian, classical solutions, and quantisation 
procedure for / are identical to the g.w. case that we recalled above, so that 
his variable / is our h for gravitational waves. His Lagrangian (2.12) for 
density perturbations coincides in structure with our Lagrangian (22), and 
we discuss one and the same observable quantity (. 

It is worthwhile to quote explicitely the attempted rigorous proof [39]: 
"Since the action (2.12) also contains a factor <p/p we also have to set its 
value to the value at horizon crossing, this factor only appears in normalizing 
the classical solution. In other words, near horizon crossing we set 




(22) 
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where / is a canonically normalized field in de-Sitter space. This produces the 
well known result...". And the author immediately writes down the square 
of the "standard inflationary result", with the square of the factor 0*/p* in 
the denominator of the final expression. 

Let us try to traverse in practice the path to the "well known result" . (To 
be fair to the author, the derivation of the "standard inflationary result" does 
not appear to be the main purpose of his paper [39], so my criticism does 
not imply anything about other statements of that paper.) The factor <p/p in 
(2.12) of the cited paper is our factor ^ in Eq.(22). It is recommended [39] 
to combine the results for the g.w. variable h with the prescription £ = -j=h. 
So, instead of Eq.(15), we would have to write 



q = C 



'% CLq 1 



k e -in(??-?7o) _|_ \J e in(v-Vo) 



2 a 

The canonically conjugate momentum seems to be 

d(' n \a J 



(24) 



(25) 



The time derivative of 7 should be neglected, as 7 is either a constant or 
a slowly changing function at times near r] . Therefore, we would have to 
write, instead of Eq.(16), the following relationship: 




_\y e -Mv-Vo) _|_ \)\ £ in(v-Vo) 



(26) 



The commutation relations are given by 

[q, p] = ih, 



b, b f 



1. 



One is encouraged and tempted to think that the quantum state |0 S ), 
annihilated by b, namely 

b|0 s ) = 0, 

is the vacuum state of the field (, i.e., the ground state of the Hamiltonian 
associated with the Lagrangian (22). The calculation of the mean square 
value of ( at r\ = i] produces the result 

(0 s |q 2 |0 s ) = --, 
2 7o 
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in which the 'zero in the denominator' factor is manifestly present and 
squared, as the "well known result" prescribes. 

In the limit of very small y/j one obtains the divergence of initial ampli- 
tudes, which is in the heart of all inflationary predictions. (In the published 
version [40] of the e-paper [39], the road to the "well known result" recom- 
mends, possibly due to a misprint, the diametrically opposite prescription 

p 

which would send the factor 7 to the numerator of the above calculation. It 
looks as though the 'rigorous' inflationary predictions fluctuate between zero 
and infinity.) 

In inflationary literature, the power spectrum P-ji(k) of curvature pertur- 
bations is usually written in the form 

Pn(k) = 

where 



and Uk are the mode-functions {u k = fi n in our notations) satisfying Eq.(19) 
with the initial conditions 

u k = -^=e~ ikn for 7] -> -00. (27) 




2k 

As one can see from the expression for P-ji(k), in inflationary theory, 
which is based on the initial conditions (27), the divergence of Pn(k) in the 
limit of very small ^7 is present from the very beginning of the evolution 
of the perturbations. To put it differently, the divergence takes place from 
the very early high-frequency regime, where by the physical meaning of the 
problem we were supposed to choose a minimal amplitude of the 'gauge- 
invariant' metric perturbation ( (or, in other words, a minimal amplitude of 
the curvature perturbation (). 

The crucial point of our discussion is that the temptation to interpret |0 S ) 
as the vacuum state for the field ( is, in fact, a grave error. The calculation 
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of the mean-square value of the canonically conjugate momentum p gives 

(0 s |p 2 |0 s ) = ^70, 
so that the factor y/j cancels out in the uncertainty relation 

A A & 

AqAp = -. 

The derived numbers clearly indicate that the quantum state |0 S ) is not 
a genuine (ordinary) vacuum state |0) for the dynamical variable (, but, on 
the contrary, is a multi-particle (squeezed vacuum) state. This is why we 
have used the subindex s. 

To show how the states |0) and |0 S ) are related, we shall first transform 
the operators. Let us introduce the annihilation and creation operators c, 
according to the Bogoliubov transformation 

b = uc + vc\ b f = u*c ] + v*c, (28) 

where 

u = cosh r, v — e t2<t> sinh r. (29) 

The parameters r and are called squeeze parameters. 
Let us assign the following values to r and 0: 

e 2r = 7 , = n(r) - r) ) or e~ 2r = 7, = n(r? - 770) + -. (30) 

We shall now use the substitution (28), together with Eqns (29) and (30), to 
Eqns (24) and (26). The factor 1/^/7 cancels out in Eq.(24) and the factor 
y/ 7 ) cancels out in Eq.(26). In terms of c, c 1 ", the operators q, p will take the 
final form: 



f h a 



ce -in(v~Vo) _|_ c t e »"(»?-»?o) 
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(31) 



(32) 



2 a 

The genuine vacuum state for the variable ( (i.e. the ground state of the 
corresponding Hamiltonian) is defined by the condition 

c|0) = 0. 
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Calculating the mean square values of q and its canonically conjugate mo- 
mentum p, we find 



(0|q 2 |0) = (0|p 2 |0) 



h 



AqAp = -, 



as it should be. 

Taking into account the definition (21), we finally derive the initial rms 



Extrapolating the initial time 770 up to the boundary between the adiabatic 
and superadiabatic regimes at r\ = rji, we arrive at the estimate 



This evaluation, plus the constancy of the quantity fx/a^/j throughout the 
long- wavelength regime, is the foundation of the result according to which the 
final (at the end of the long-wavelength regime) amplitudes of gravitational 
waves and density perturbations should be roughly equal to each other [33]. 

There is no dimensional parameter which could be regulated in such a 
way as to make one of the amplitudes several orders of magnitude larger 
than another. In terms of the 'scalar' and 'tensor' metric amplitudes this 
means that hx/hs ~ 1 for all 7's. More accurate calculations along the same 
lines produce C = \f2Anlpi in expression (7) for density perturbations. 

Certainly, the correct quantisation procedure (31), (32), as opposed to 
the incorrect (inflationary) procedure (24), (26), could be formulated from 
the outset of quantisation. Mathematically, the Lagrangians (13) and (22) 
are alike, if in (13) one means a by a, and replaces h with (. 

The derivation of the Hamiltonian for S-perturbations repeats the deriva- 
tion for gravitational waves. Using the canonical pair q — p — dL/d/i' for 
/is, w e arrive at the Hamiltonian (compare with Eq.(98) in Ref. [33]) 



where the coupling to the external field is now given by the function 17(77) = 



value of the variable ( = n/a^fy: 




(33) 




H(t7) = ne'e + ctc 1 " + a*c 2 , 



(34) 



(i/2)(5'/o). 
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The Heisenberg equations of motion for the Heisenberg operators c(r/), 
c^(t/) lead to classical equation (19). The asymptotic expressions for the 
Heisenberg operators, 

c(n) = ce~ in ^- m \ c\rf) = c^e in( ^°\ 

are participating in Eqs. (31), (32). Clearly, the vacuum state |0), defined 
as c(r])\0) = 0, minimises the oscillator's energy (34). 

Since at times near t)q the coefficients a/a^ and a/a® are close to 1, the 
equality of the initial values for h rms and ( rms follows already from the simple 
comparison of the Lagrangians (13) and (22). 

The relationship between the above-mentioned genuine vacuum state |0) 
and the squeezed vacuum state |0 S ) is determined by the action of the squeeze 
operator S(r,(f>) on |0): 

\0 8 )=S{r,<j>)\0), 



where 

S(r, 0) = exp 



1 

2' 

The mean number of quanta in the squeezed vacuum state is given by 



(0 s |c T c|0 s ) = sinhV = 



2_ 1-7 



2^7 



This is a huge and divergent number, when the 'zero in the denominator' 
factor goes to zero. Therefore, the "standard inflationary result" for S- 
perturbations is based on the wrong initial conditions, according to which 
the initial amplitude of the (^-perturbations can be arbitrarily large from the 
very beginning of their evolution. 

Moreover, the initial amplitude is assumed to go to infinity in the most 
interesting limit of v /7 — > and n — > 1. If ^7 does not deviate from 1 
too much, then the mean number of quanta in the squeezed vacuum state 
is acceptably small, and the wrong initial conditions give results sufficiently 
close to the correct ones. However, as in the Landau example mentioned 
above, if the wrong formula gives acceptable answers for some range of x, 
this does not make the wrong theory a correct one. (Finally, if y/j = 1, then 
a(t) oc t, a{rj) oc e v , w = —1/3. From this model of cosmological evolution 
the study of relic gravitational waves has began in the first paper of Ref.[2].) 
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In terms of the classical mode functions, it is the function that 
should satisfy the classical version of the initial conditions (31), and not the 
function fi/a, which is postulated by the inflationary requirement (27). They 
both are so-called 'gauge-invariant' variables, but their physical meaning is 
drastically different. The original derivations of the "well-known result" were 
guided simply by the visual analogy between the function u — n in theory of 
density perturbations and the function \i in the theory of gravitational waves 
already developed by that time. 

The assumption of arbitrarily large initial amplitudes of curvature per- 
turbations or, technically speaking, the choice of the initial multi-particle 
squeezed vacuum state |0 S ) for (, instead of the ordinary vacuum sate |0), is 
the origin of the absurd "standard inflationary result" . Certainly, this wrong 
assumption cannot be the basis of observational predictions for cosmology. 

5 Conclusions 

The grossly incorrect predictions of inflationary theorists should not be the 
reason for doubts about the existence and expected amount of relic gravi- 
tational waves. The generation of relic gravitational waves is based on the 
validity of general relativity and quantum mechanics in a safe cosmologi- 
cal regime where quantisation of the background gravitational field is not 
necessary. 

In our numerical evaluations, we also assumed that the observed large- 
angular-scale anisotropies of CMB are caused by cosmological perturbations 
of quantum-mechanical origin. This is not necessarily true, but it would be 
quite disastrous if it proved to be untrue. 

It is quite a challenge to imagine that the natural and unavoidable quantum- 
mechanical generation of cosmological perturbations is less effective than 
anything else. In any case, if relic gravitational waves are not discovered at 
the (relatively high) level described in this contribution, the implications will 
be much more serious than the rejection of one inflationary model or another. 
The reality of our time is such that if the proposal is not properly 'sexed-up', 
it is not very likely to be funded. But the ultimate truth lies in the fact that 
the real physics of the very early Universe is much more exciting than the 
artificial hullaballoo over popular words such as 'inflation' or 'inflationary 
gravitational waves'. 
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Hopefully, relic gravitational waves will be discovered in experiments, 
which are already in the well-developed stage. I personally would think that 
this is likely to happen first in dedicated ground-based observations, such 
as the recently approved Cardiff-Cambridge-Oxford collaboration CLOVER 
[41]. Let us hope this will indeed be the case. 
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